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Abstract 


33 

34 The Air Force Weather Agency (AFWA) - NASA blended snow-cover product, called 

35 ANSA, utilizes Earth Observing System standard snow products from the Moderate- 

36 Resolution Imaging Spectroradiometer (MODIS) and the Advanced Microwave 

37 Scanning Radiometer for EOS (AMSR-E) to map daily snow cover and snow-water 

38 equivalent (SWE) globally. We have compared ANSA-derived SWE with SWE values 

39 calculated from snow depths reported at -1500 National Climatic Data Center (NCDC) 
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40 co-op stations in the Lower Great Lakes Basin. Compared to station data, the ANSA 

41 significantly underestimates SWE in densely-forested areas. We use two methods to 

42 remove some of the bias observed in forested areas to reduce the root-mean-square 

43 error (RMSE) between the ANSA- and station-derived SWE. First, we calculated a 5- 

44 year mean ANSA-derived SWE for the winters of 2005-06 through 2009-1 0, and 

45 developed a five-year mean bias-corrected SWE map for each month. For most of the 

46 months studied during the five-year period, the 5-year bias correction improved the 

47 agreement between the ANSA-derived and station-derived SWE. Flowever, anomalous 

48 months such as when there was very little snow on the ground compared to the 5-year 

49 mean, or months in which the snow was much greater than the 5-year mean, showed 

50 poorer results (as expected). We also used a 7-day running mean (7DRM) bias 

51 correction method using days just prior to the day in question to correct the ANSA data. 

52 This method was more effective in reducing the RMSE between the ANSA- and co-op- 

53 derived SWE values, and in capturing the effects of anomalous snow conditions. 

54 

55 Introduction 

56 

57 Significant reductions in the extent of spring snow cover in the Northern Hemisphere 

58 have been measured during the satellite era (Brown, 2000; Brown et al., 2010) and over 

59 the past 90 years, including about 50 years of the pre-satellite era (Brown and 

60 Robinson, 201 1 ). Accelerated warming is most likely responsible for the observed 

61 accelerated decrease in Northern Hemisphere spring snow cover in recent years 
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62 (Brown and Robinson, 2011). It is important to monitor snow depth and snow-water 

63 equivalent (SWE) as well as snow extent. 

64 

65 Consistent and reliable snow maps are needed by the climate-modeling community to 

66 to improve the predictive capabilities of the models, and to validate their performance. 

67 Several global snow maps are in widespread use (e.g., Robinson et al., 1993; Ramsay, 

68 1998; Hall et al., 2002; Kelly et al., 2004; Frei et al., submitted). Most satellite-derived 

69 global snow maps have relied on utilization of data from one primary satellite, though 

70 the NOAA National Ice Center utilizes multiple satellite datasets as well as station data 

71 to construct daily maps (Helfrich et al., 2007). 

72 

73 The Great Lakes Basin is a challenging study area for mapping snow using satellite 

74 data for several reasons. Persistent cloud cover is a serious issue when observing 

75 snow using sensors operating in the visible-near-infrared (VNIR) wavelengths. And 

76 passive-microwave algorithms have limitations when mapping the shallow often 

77 ephemeral snow in this region as well as the wet snow - daytime temperatures are 

78 typically above 0°C in the lower Great Lakes area. Moreover, the large footprint size of 

79 passive-microwave sensors means that small scale lake-effect snow events may be 

80 missed. 

81 

82 A blended snow-cover product has been developed jointly by the U.S. Air Force 

83 Weather Agency (AFWA) and the Hydrospheric and Biospheric Sciences Laboratory at 

84 NASA / Goddard Space Flight Center. A detailed description of the product, derived 
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85 from the AFWA - NASA Snow Algorithm (ANSA) may be found in Foster et al. (201 1 ). 

86 The objective of the present work is to characterize the accuracy of the SWE derived 

87 from the ANSA snow maps in a region within the Lower Great Lakes Basin for the 

88 winters of 2005-06 through 2009-10 using meteorological-station data. We also 

89 describe methods to improve the accuracy of the ANSA-derived SWE by reducing the 

90 bias in the ANSA measurements. 

91 

92 

93 Study area 

94 

95 The Great Lakes system (Figure 1 ) is the largest surface freshwater system on the 

96 planet. It drains an area of approximately 1 ,600,000 km 2 and stretches from northern 

97 Minnesota / western Ontario to the Gulf of St. Lawrence. The mean annual flow, as 

98 measured at the mouth of the St. Lawrence River, is approximately 12,600 m 3 /s/km. 

99 Snow cover is prevalent during winter, and snowfall averages more than 80 cm in 

100 southern locations and more than 250 cm at a few locations in the lee of Lake Superior 

101 and Lake Ontario and at the highest elevations. Basin relief is rather low - on the order 

102 of hundreds of meters. Vegetation consists of transitional mixed forests, northern 

103 hardwoods, and Great Lakes spruce and pine forest. 


104 
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108 The Great Lakes Basin is subject to lake-effect snowfall, particularly elevated areas in 

109 the lee of the Great Lakes. Lake-effect snowfall is produced when cold winds move 

110 across long expanses of relatively-warmer lake water; the warmer water provides 

111 energy and a source of water vapor. Snow deposited on the leeward shores can 

112 accumulate to significant depths in relatively short periods -- 0.5 m or more in 24 hours. 

113 Though the heaviest accumulations typically occur within 80 km of the lakes, on 

114 occasion lake-effect snow may fall 320 km downstream. For example, the highlands of 

115 West Virginia receive most non-storm snowfall from lake-effect snows from Lake Erie. 

116 

117 

118 Air Force - NASA Snow Algorithm (ANSA) snow product 

119 

120 An example of the ANSA snow product is provided in Figure 2. The product utilizes the 

121 Moderate-Resolution Imaging Spectroradiometer (MODIS) standard snow-cover maps 

122 (Hall and Riggs, 2007) and the Advanced Microwave Scanning Radiometer for EOS 

123 (AMSR-E), a passive-microwave instrument, standard snow-water equivalent (SWE) 

124 maps (Kelly et al., 2004; Kelly 2009: Tedesco, 201 1 ) to map daily snow cover and SWE 

125 globally. These products have been described in great detail elsewhere therefore only 

126 a brief description will be provided here. 


127 



7 


128 

129 

130 

131 

132 




Snow by both sensors 

Snow by AMSR E, MODIS cloud or no data 

Snow by MODIS, AMSRE no snow or orbit gap 

No snow by MODIS or AMSR E but cloud obscured 

No snow: no snow by MODIS in clear view but, AMSR E detects snow 

Cloud by MODIS in AMSR E orbit gap 

Snow tree land by both MODIS and AMSR E 


Figure 2 


MODIS standard snow maps (MOD10C1) provide high-quality, daily and global 
cover maps at a spatial resolution of up to 500 m (Hall et al., 2002; Riggs et a!., 
A 500-m resolution MODIS snow-cover map is shown in Figure 1b, however, in 


snow- 

2006). 

the 


133 


ANSA product, we use 5-km resolution snow products as the default for mapping snow 
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134 cover. AMSR-E data can map snow through clouds and darkness and provide 

135 estimates of SWE at a spatial resolution of ~25 km. AMSR-E-derived snow cover is 

136 used when clouds preclude MODIS from providing a snow map. 

137 

138 SWE derived from the ANSA product comes only from the AMSR-E. The MODIS VNIR 

139 bands cannot directly measure snow depth, nor can they image through cloud cover 

140 which is persistent in the Lower Great Lakes region during the wintertime. 

141 

142 Previous work (Figure 3) has shown that use of the ANSA product enables improved 

143 mapping of snow-cover extent in the Lower Great Lakes region relative to using either 

144 MODIS or AMSR-E maps alone (Hall et al. , 2009). Use of the ANSA snow products 

145 was also found to improve the mapping of snow-cover extent for the 2007-08 winter in a 

146 mountainous area in the eastern part of Turkey where the elevation ranges between 

147 850 and 3000 m (Akyurek et al., 201 0). 91 % agreement was obtained between the 

148 ANSA snow maps and in-situ observations for February 2008. This is the first time that 

149 the ANSA snow cover product was evaluated in a mountainous area. Daily snow data 

150 collected at 36 meteorological stations were used for validation. 

151 
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Figure 3 


Methodology 

We use National Climatic Data Center (NCDC) co-op station data (Figure 4) for the five 
winters from 2005-06 through 2009-10 in the Lower Great Lakes region, to compare 
with ANSA-derived SWE (Figure 4). Co-op snow depth data were interpolated to 
develop a daily map (see a sample of an interpolated map for 1 December 2007, in 
Figure 5) and then converted to SWE using two different density conversion factors: 0.2 
and 0.3 (representing snow densities of 0.2 and 0.3 g/cm 3 , respectively). Since the 
exact snow density was unknown (snow density is not routinely made at co-op stations), 
initially, we converted snow depth to SWE using the above snow densities, which are 
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166 reasonable for the conditions and time of year. (In reality the snow density changes 

167 over time as the snow metamorphoses.) 

168 
169 



170 


Figure 4 
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No snow Deeper snow 

Figure 5 


Daily difference maps were then constructed to evaluate the accuracy of ANSA-derived 
SWE as compared with SWE derived from interpolated station data. These difference 
maps (ANSA SWE minus station-derived SWE) can be considered a measure of 
deviation (or error) of the ANSA SWE from the “truth,” the co-op station data. The 
RMSE quantifies an over- or under-estimate of the actual SWE. RMSE was calculated 
for each pixel for each day’s snow map. RMSE was averaged to get a daily value for 
the entire region of interest (or domain). 
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181 We then developed two methods to bias correct the ANSA daily SWE maps. The 

182 attributes of each method are described below. 

183 

184 First, a monthly error was determined by calculating the mean difference of ANSA SWE 

185 minus station-derived SWE for each cell for each day of each month (November, 

186 December, January, February and March of each year). This calculation was repeated 

187 for each of the five winters, 2005-06 through 2009-1 0, resulting in a 5-year bias map for 

188 each month. To remove the 5-year bias from each daily ANSA map, the mean-monthly 

189 difference was subtracted from the ANSA SWE in each cell; this produced bias- 

190 corrected SWE maps for each day of each month studied. 

191 

192 Second, a 7-day running mean (7DRM) using days immediately before the day of 

193 interest was calculated and the average difference map of those seven days was used 

194 as a bias to correct the ANSA SWE data on the following day. For example, the 7-day 

195 period, 1 -7 January, is used to calculate the 7DRM to bias correct the ANSA SWE on 

196 the 8 January snow map, then a correction would be calculated for 2-8 January, from 

197 which a new bias is used to correct the 9 January map, and so on. The daily RMSE, 

198 the difference between the ANSA data corrected with the 7DRM method, and the 

199 corresponding station data, is computed each day. 

200 
201 

202 Results & Discussion 


203 
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The difference between the SWE derived from the interpolated co-op data compared to 
the ANSA-derived SWE for each cell without any bias correction is shown in Figures 6a 
and 6b. The blue colors indicate where the ANSA map underestimates SWE as 
compared to station data, and the pink colors indicate that ANSA overestimates SWE. 
The areas of best agreement are shown in neutral colors. 
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Figure 6 
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212 

213 In January and February of 2008, ANSA underestimates SWE in densely-forested areas 

214 such as in the Upper Peninsula of Michigan and in the Adirondacks in New York State 

215 by up to ~75 mm (see red circles in Figure 6a). Even when using its forest-fraction 

216 adjustments (Foster et al., 2005), passive-microwave algorithms still underestimate 

217 SWE in dense forests. 

218 

219 Conversion of snow density to SWE. What snow density should be used to convert 

220 the co-op snow depths to SWE? Ideally the snow density should vary as described in 

221 Foster et al. (2005). Flowever, lacking specific information on snow-cover 

222 metamorphism in this study area, we decided to use a fixed snow density to compute 

223 the snow depth from the AMSR-E SWE measurements. In both the January and 

224 February 2008 difference maps (Figures 6a & b), there is overall better agreement when 

225 a snow density of 0.3 g / cm 3 was used to convert the snow depths to SWE versus using 

226 0.2 g / cm 3 (Figures 6c & d). Note the substantial differences when different snow 

227 densities are used to convert snow depth to SWE. This demonstrates the need to 

228 improve density estimates, especially in areas where snow conditions change rapidly, 

229 and to incorporate dynamic features into SWE algorithm so that the density changes 

230 with time. 

231 

232 

233 Improvement of ANSA SWE calculation using bias correction. There are both 

234 systematic (bias) and random errors associated with the passive-microwave 
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235 measurements. In order for the remotely sensed SWE observations to be useful for 

236 climate modelers, for instance, it is necessary to have both an unbiased SWE estimate 

237 and a quantitative, rather than qualitative, estimate of the uncertainty (random errors). 

238 This is a critical requirement for successful assimilation of snow observations into land 

239 surface models. 

240 

241 First, we experimented with the 5-year bias information to produce “5-year bias- 

242 corrected” SWE maps as shown in Figure 7a. Only areas that are currently snow 

243 covered are displayed on the maps, and non-snow-covered terrain is shown in green. 

244 Using the 5-year bias correction, the RMSE improves on most days in most months. 

245 However when anomalous snow conditions occur, the 5-year bias correction can 

246 substantially increase the RMSE. 

247 

248 In general, there is considerable improvement in the accuracy of the ANSA SWE 

249 measurement as compared to station data, when we remove the 5-yr average bias. 

250 This can be seen by comparing Figures 7a and 7b, and by noting that the RMS errors 

251 are lower when the 5-yr bias is removed as seen in Figure 8. This bias-correction 

252 technique works quite well when snow conditions on a given day are reasonably close 

253 to “average.” However when snow conditions are anomalous, as seen in Figure 9, use 

254 of the 5-yr bias correction can result in higher RMSE values (see green line). 

255 

256 


257 
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Figure 7 
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263 Between 6 and 8 February 2006 we see a change in the RMSE values (Figure 9) 

264 because the original ANSA data (with no bias correction) matches the interpolated snow 

265 depth map better than does the 5-year bias-corrected data. For most of February 2006, 

266 snow conditions differed considerably from the 5-year average and therefore the 5-year 

267 bias correction did not offer improvement in this case. For example, in Indianapolis, 

268 Indiana, which is representative of the Lower Great Lakes region, the temperatures 

269 were warmer than normal in February and the snowfall was much lower than average. 

270 Only three days recorded 2.5 cm of snow on the ground (4-6 February 2006), and 

271 there was no day that recorded more than 5 cm of snow. In addition, there was no snow 

272 on the ground (other than a trace) after 6 February. 

273 In all cases that we studied during the five winters, the 7DRM bias-correction technique 

274 worked better than using either no bias correction for the ANSA SWE data, or using the 

275 5-year average bias-correction technique. The 7DRM approach captures the snow 

276 conditions just before the day in question. Future work calls for trying additional ways to 

277 bias-correct the ANSA data using data just prior to the day in question, but allowing a 

278 few days’ delay to acquire the data to perform the bias correction. 

279 Models that improve the evolution of snowpack parameters, including grain size 

280 information, and use of dynamic algorithms that better account for changes in snow 

281 density, should be the focus of future work. A dynamic method of calculating snow 

282 density in a passive-microwave SWE algorithm would likely result in improved SWE 

283 estimates from microwave sensors. 


284 
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285 

286 Conclusions 

287 

288 We have examined the ability of the ANSA blended-snow product to measure SWE in 

289 the Lower Great Lakes region of the U.S. The ANSA product underestimates SWE in 

290 dense forests. This is a known limitation of passive-microwave SWE algorithms. Thus 

291 we employed bias-correction methods to reduce the errors when measuring SWE from 

292 ANSA. First, we calculated the 5-year mean-monthly bias (difference between co-op- 

293 derived SWE and AMSR-E SWE) from the daily maps from each month (November 

294 through March) of each winter from 2005-06 to 2009-10. Those values, on a cell-by-cell 

295 basis, were subtracted from the ANSA SWE values for each day of each corresponding 

296 month. Results show an improvement in agreement between the co-op station-derived 

297 and ANSA SWE after the monthly bias was removed. 

298 

299 Estimation of snow density is needed to convert the co-op station snow depths to SWE 

300 values, so that ANSA SWE can be compared with “ground truth.” For our study area, 

301 use of a conversion factor of 0.3 (corresponding to a snow density = 0.3 g/cm 3 ) provides 

302 overall better agreement between ANSA and co-op SWE for both January and February 

303 2008, than when we used a snow density of 0.2 g/cm 3 . 

304 

305 For most of the months during the five-year period, the 5-year bias correction improved 

306 the agreement between the ANSA-derived and station-derived SWE. However, 

307 anomalous months such as when there was very little snow on the ground compared to 
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308 the 5-year mean, or months in which the amount of snow was much greater than the 5- 

309 year mean, showed poorer results. We also used the 7DRM bias-correction method 

310 using the 7 days just prior to the day in question. As before, we then corrected the 

311 ANSA data. This method was more effective in reducing the RMSE between the ANSA- 

312 and co-op-derived SWE values, and in capturing the effects of anomalous snow 

313 conditions. 

314 
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390 Figure captions 

391 

392 Figure 1 a. Terra MODIS image acquired on 1 6 February 2008 showing snow cover in 

393 the vicinity of the Great Lakes in the northern United States and southern Canada. 

394 There is ice cover on Lake Erie but the other Great Lakes are mainly cloud covered. 

395 Image courtesy of the Earth Observatory Image of the Day 

396 http://earthobservatory.nasa. qov/IOTD/view.php?id=8485 . Figure 1b. MODIS snow- 

397 cover fraction (SCF) map of the same area shown in Figure 1 a. 

398 

399 Figure 2. ANSA blended-snow product for 26 January 2007 in Lambert Azimuthal polar 

400 projection (Foster et al., 201 1 ). 

401 

402 Figure 3. Relationship of the Percent of Agreement of the ANSA product, and the 

403 MODIS and AMSR-E input products, alone, as compared to meteorological station data 

404 for the lower Great Lakes region for mapping snow-cover extent in 2003. The AMSR-E 

405 contribution (green) becomes more important in late February during periods of 

406 cloudiness when MODIS cannot map the snow (after Flail et al., 2009). 

407 

408 Figure 4. Dots represent locations of NCDC co-op stations used in this study. 

409 

410 Figure 5. 1 Dec. 2007 snow depth map interpolated from co-op station measurements. 

411 Snow depths are shown in various shades of black and white (lighter grey indicates 

412 deeper snow). 
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413 

414 Figure 6. Mean difference between ANSA and station-derived SWE for January (A & B) 

415 and February 2008 (C & D) using 0.2 conversion factor (A & C), and 0.3 conversion 

416 factor (B & D). Conversion factors are used to convert co-op station snow-depths to 

417 SWE and are based on snow density. The blue colors indicate where the ANSA map 

418 underestimates SWE as compared to station data, and the pink colors indicate that 

419 ANSA overestimates SWE. 

420 

421 Figure 7a. Mean difference between ANSA and station-derived SWE for 2 January 

422 2009. Figure 7b. Mean difference between ANSA and station-derived SWE for 2 

423 January 2009, but in this image, the 5-year bias was subtracted from the result in Fig. 

424 7a, for each cell. In these images, the non-snow-covered land areas are shown in 

425 green. 

426 

427 Figure 8. Plot of the RMSE for January 2009. The top (blue) line corresponds to Figure 

428 7a (no bias correction), and the green line corresponds to Figure 7b (5-year bias 

429 correction method). The orange line represents results using the 7-day running mean 

430 (7DRM) bias-correction method. 

431 

432 Figure 9. Plots of the RMSE for the month of February 2006. The blue line represents 

433 the case where the original ANSA data were compared with the station-derived SWE 

434 (no bias correction), and the green line represents the case where the 5-yr bias 

435 correction was used. Note in the middle of the month that the errors are greater when 
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the 5-yr bias correction is used. The orange line represents results using the 7-day 
running mean (7DRM) bias-correction method. 
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